Mimicking a turbulent signal: sequential multiafRne processes 
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An efficient method for tlie construction of a multiaffine process, witfi prescribed scaling expo- 
nents, is presented. At variance witfi tlie previous proposals, this method is sequential and therefore 
it is the natural candidate in numerical computations involving synthetic turbulence. The applica- 
tion to the realization of a realistic turbulent-like signal is discussed in detail. The method represents 
a first step towards the realization of a realistic spatio-temporal turbulent field. 
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In recent years the relevance of multifractal measures 
and multiaffine processes in many fields (mainly fully de- 
veloped turbulence) has been well understood 0-^- In 
different contexts, for instance numerical simulations and 
comparison of theoretical models with experimental data, 
a rather natural problem is the construction of artifi- 
cial signals mimicking real phenomena (e.g. turbulence). 
In particular it is important to have efficient numerical 
techniques for the construction of a multiaffine field (^{x) 
whose structure functions scale as 



{\4>{x + r) - 4>{x)\'') 



(1) 



where (• • •) indicates a spatial (or temporal) average, r 
varies in an appropriate scaling range and the exponents 
C,q are given. The most interesting case, and the most 
physically relevant, is when C,q is a nonlinear function of 
q, that is a strictly multiaffine field. 

Let us first notice that the generation of a multiaffine 
function is much more difficult that the generation of 
a multifractal measure, which can be obtained with a 
simple multiplicative process generalizing the two scales 
Cantor set. 

Up to now, there exist well established methods for the 
construction of multiaffine fields |^-|| , see for a short 
review. All of these methods share the common charac- 
teristic of being not sequential: the process is build as a 
whole in an interval (in space or time) of fixed length. To 
extend the interval one has to rebuild the process from 
the beginning. This is an evident limitation if one is 
interested in constructing a temporal signal mimicking, 



for example, those obtained by an anemometer measure- 
ment. Furthermore, non-sequential algorithms require 
always a huge amount of stored data. 

In this letter we introduce a simple and efficient se- 
quential method for the construction of a multiaffine 
function of time u{t) with prescribed statistical proper- 
ties. The guideline of our approach will be the repro- 
duction of a turbulent-like temporal signal. Though the 
basic idea on the construction of the multiaffine process 
comes from fully developed turbulence, nevertheless the 
method is general and can be applied to any signal. 

A typical anemometer measurement gives a I- 
dimcnsional string of data representing the one-point tur- 
bulent velocity u{t) along the direction of the mean flow 
U . According to the Taylor hypothesis Q, for small tur- 
bulence intensities u <^ U, the time variations of u can 
be assumed to be due to the advection (with velocity U) 
of a frozen turbulent field past the measurement point, 
so that 

Su{t) — u{x, t + t) — u{x, t) = 

= u{x -UT,t) ~u{x,t) = 5u{l) (2) 

where £ — Ut. Therefore, once the spatial scaling (|l|) is 
given, we have: 



Sq{r) = {\u{t + r)-u{t)n^r':^ 



(3) 



The frozen field is the result of the superposition of tur- 
bulent patterns (eddies) of many different sizes i, whose 
contribution to the time variation of the velocity decays 
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with a typical correlation time Tsweep ~ £/U. For the 
sake of simplicity, in the following, we shall introduce a 
set of reference scales = 2^" at which scaling proper- 
ties will be tested. With this picture in mind, we repre- 
sent the signal u{t) by a superposition of functions with 
different characteristic times, representing eddies of var- 
ious sizes 



N 



N 



(4) 



The functions Vn{t) are defined via a multiplicative pro- 
cess 



Vn{t) = gn{t)xi{t)x2{t) . . . Xn{t) , 



(5) 



where the gn{t) are independent stationary random 
processes, whose correlation times are the sweeping 
timescales t„ = £n/U — 2~" (assuming U — 1) and 
(5n) = where h is the scaling exponent. For fully de- 
veloped turbulence h = 1/3. Scaling will show up for all 
time delay larger than the UV cutoff tn and smaller than 
the IR cutoff Ti. The Xj{t) are independent, positive de- 
fined, identical distributed random processes whose time 
correlation decays with characteristic time Tj . The prob- 
ability distribution of xj determines the intermittency of 
the process. 

The origin of (||) is fairly clear in the context of fully de- 
veloped turbulence. Indeed according to the Refined Sim- 
ilarity Hypothesis of Kolmogorov , we can identify 
Vn with the velocity difference at scale £n and Xj with 
(ej/ej-iY^^, where ej is the energy dissipation at scale 

It is easy to show, with a simple argument, that 
the process constructed according to (^,^) is multiafhne. 
Because of the fast decrease of the correlation times 
Tj = 2~^, the characteristic time of Vn{t) is of the or- 
der of the shortest one, i.e., t„ — 2^". Therefore, the 
leading contribution to the structure function Sq (r) with 
T ^ Tn will stem from the n-th term in (jj). This can 
be understood nothing that in the sum u{t + t) — u{t) ~ 
+ r) — Vk{t)] the terms with k < n are negli- 
gible because Vk{t + T) ~ Vk(t) and the terms with k > n 
are subleading. Thus one has: 

SM ~ {\vn\'} ^ (|g„|«)(x«)" ^ r^-i°«^<-'> (6) 

and therefore for the scaling exponents (^) 

Cq = hq - \og^{x^) . (7) 

The limit of an affinc function can be obtained when all 
the Xj are equal to 1. 

The above results can be proved in a more rigorous way 
considering the second order structure function 5*2 (t). 
Using the definitions (^,^ and stochastic independence 
one obtains: 



52(t) = 2j2[Mtf) - {Vnit)Vnit + T))]. (8) 
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Let us now introduce the normalized correlation func- 
tions for gnit) and Xj{t) 



C 



F 



s \ _ {xj{t + s)xj{t)) 



(9) 



(10) 



Plugging into (^ the definition (H) one obtains 



N 



Tl = l 



52(r) = 2Y,{gl){xY 1 - C(-)F(-) • • • Fi- 

^ V Tn Ti Tn 



(11) 



By shifting the summation index in the above expression, 
n — > n — 1, one obtains for r 1, 



S2{2T)^2''^{xY'S2{r) 
which leads to the scaling behavior 



(12) 



52 (r) 



with C2^2h- log2(a;^) . (13) 



A similar computation can be performed for the higher 
order structure functions. The generic Sq{T) can be 
expressed as a linear combination of terms scaling as 
T-Cmi . . . ^Cmfc -vjfith mi + . . . + mfe = q. From the con- 
vexity of Cg it follows that the leading contribution 
to (r) for small t is given by Sq{T) ^ t^" , with the 
exponents C,q as defined in (^. 

The key point in the above arguments is that the dom- 
inant contribution to the structure function Sq{T) comes 
from octaves n such that t„ ^ r, that is locality. 

The constraints for locality can be captured with a sim- 
ple argument. Indeed for t„ <C t we have that {\vn(t + 
T)—Vn{t)\'^) ^ {\vn\'^) ^ 2^"'''', therefore UV Convergence 
requires (^q > 0. Similarly, when r„ ^ r we have that: 

{\v„{t + t) - Vnim ~ ir/Tn)''/'{\Vn\'^) ^ 2-<^'^-ll^\ 

for stochastic processes such that C{x) = 1 — 0(x) and 
F(x) = 1 — 0[x). Therefore, convergence in the lat- 
ter case requires < 9/2. We observe that the last 
condition is different from the usual locality condition 
Cq < q [|l3| which holds for differentiable processes where 
C{x) = 1 - 0{x^) and F{x) = 1 - 0{x^). 

Regular behavior for very short time delays Su{t) ^ t, 
physically related to the presence of dissipation, can be 
simply achieved in our model by smoothing gnit) and 
Xn{t) over a time interval smaller then the UV cutoff t^. 

The numerical implementation of the method proposed 
above is very simple. The stochastic process Xj{t) can be 
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easily generated via the nonlinear Langevin differential 
equations: 



I dV ^ 2 
axj — — at + . — aWj 

Tj dXj V Tj 



(14) 



where V{x) — oo for x < a {a positive constant) and 
V{x) — > oo for X — > oo. It is clear that the Xj so obtained 
have the same probability density function independently 

of Tj. ^ 

Similarly for the gn one can use the evolution law 



1 dY , 
dgn = 3 — dt 

Tn dgn 



dWn 



(15) 



where Y{g) — > oo as |g| — > cxd and (T„ = i^- 

Numerical tests have been performed adopting for the 
stochastic differential equations (|l^,|l5|) the following po- 
tentials: 



V{x) 



-21nx 



for (1 - 6)1/3 <x <{1 + 6)1/^ (16) 



1/3 



and V{x) — oo otherwise, where < 6 < 1 , while 



Y{g) 



1 



-,9 



(17) 



implemented by using any hierarchical constructions pre- 
viously proposed |8|. Nevertheless, this way to glue 
together spatial and temporal multiafhne fluctuations 
would not be realistic, due to the absence of a real sweep- 
ing of small scales by large scales. This is connected 
to the fact that in our temporal signal, the Taylor hy- 
pothesis is introduced by hands without any real direct 
dynamical (stochastic) coupling between large and small 
scales. 

These difficulties in reproducing an Eulerian spatio- 
temporal field are absent if one considers the velocity 
statistics in quasi-lagrangian coordinates ||l^. In this 
framework a pure temporal signal would correspond to 
the velocity field felt in the moving reference frame at- 
tached to a fluid particle. The sweeping effect is thus 
removed and the characteristic time scales are the dy- 
namical eddy turnover times. Work in this direction is 
in progress. 

We thank D. Picrotti for useful discussions in the 
early stage of the work. This work has been partially 
supported by the INFM (Progetto di Ricerca Avanzata 
TURBO). 



For h — 1/3, this choice insures that Ca = 1 according to 
the scaling prescribed by Kolmogorov's law. The param- 
eter b tunes the intermittency of the signal: when b I 
we recover an affine process. 

In figure 1 we show the the quantity vf^it) which can 
be considered as the energy density dissipation of the 
turbulent signal. As one can see high intermittency is 
detected. 

The theoretical and numerical scaling laws are com- 
pared in figure 2. The computed scaling exponents are 
in perfect agreement with those given by equation (|^). 
Figure 3 shows the probability density function of the 
velocity differences Su{t) — u{t + t) — u{t) for different 
T. At large r 1 the pdf is nearly Gaussian, whereas at 
small delays the pdf is increasingly peaked around zero 
with high tails corresponding to large fluctuations with 
respect to their rms value. If one wants the process u(t) 
to have a nonzero skewness, as in turbulence, Y{g) must 
be chosen as an asymmetric function, see |^ for a suitable 
choice according to experimental data. 

In this letter we have introduced an efficient sequen- 
tial algorithm for the generation of multiaffine processes. 
This method, at variance with previous proposals, is not 
based on hierarchical construction, and can be applied to 
any multiaffine signals with specified scaling laws. Fur- 
thermore, no huge amount of memory is required for the 
numerical implementation. 

A possible, relevant, application of such a signal would 
be to use it for describing the temporal part of a syn- 
thetic turbulent velocity field. The spatial part can be 
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FIG. 1. Time series v%{t) normalized to the average for the 
model with A'^ = 15 octaves and b = 0.9. 




FIG. 2. Numerical (dots) and theoretical (line) structure 
functions 5'g(r) for the model with A'^ = 20 octaves and 
b = 0.9. The exponents are Ci = 0.39, (3 = 1,^6 = 1.65. 
The structure functions are shifted by a multiplicative factor 
for plotting purposes. 




6u / a 

FIG. 3. Probability density functions for the normalized 
velocity differences 6u{T)/a, where a = {5u^y^^, for differents 
T. For large r = 10 (b) the pdf is nearly Gaussian (dashed 
curve). For very small r = 0.001 (a) large tails are evident. 
The parameters are = 15 octaves and b = 0.9. 



4 



